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WHAT IS CLAIMED IS: 



1 . A method for locating an object relative to an array of magnetic 
sensors in an environment in which there is present a noise signal 
5 having a fundamental frequency f NO isE> the method comprising: 

generating one or more magnetic signals by means of one or 
more magnetic emitters mounted at known locations on the object, 
the one or more magnetic signals having one or more frequencies; 
during an integration time T, detecting the one or more 
10 magnetic signals and the noise signal at six or more magnetic 

detectors; and, 

determining relative amplitudes of the magnetic signals; 
wherein the one or more frequencies of the magnetic signals are 
substantially equal to frequencies at which a power spectrum of 
1 5 the detected noise signal has zeros. 



2. A method according to claim 1 wherein the integration time J is 
equal to Njlf NO!SE , where 7V 7 is an integer and Nj>2. 



20 3. A method according to claim 2 wherein the one or more 

frequencies of the magnetic signals are substantially equal to N 2 /T 
where N 2 is an integer which is not an integer multiple of N 2 . 



4. 

25 



A method according to claim 1 wherein detecting the magnetic 
signal comprises: 
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sampling outputs of the six or more magnetic detectors to 
provide corresponding sequences of samples and, for each of the 
one or more magnetic detectors; 

multiplying each of the samples by a function which varies 
periodically at a frequency at which the power spectrum of the 
detected noise signal has a zero; and, 

summing results of the multiplying. 

A method according to claim 4 wherein the function is a cosine 
function. 



A method according to claim 5 wherein detecting each of the one 
or more magnetic signals comprises, for each of the one or more 
magnetic detectors, computing the complex value A given by: 



A = 



X cos(2*^ / f s ) - S : sin(2**/; / f s ) 



or a mathematical equivalent thereof where S k m is the k^ sample 
from an m ih one of the magnetic detectors; k is an index indicating 
the order of the samples in the sequence of samples from each of 
the magnetic detectors;/^ is a frequency of the magnetic signal; f s 
is a frequency at which the samples are taken; and i is the square 
root of - 1 . 



A method according to claim 6 wherein the one or more magnetic 
signals have frequencies f q which deviate from frequencies f qi at 
which the power spectrum of the detected noise signal has zeroes 
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by a deviation amount and the method comprises performing a 
correction to the complex value A. 

A method according to claim 1 wherein the one or more magnetic 
emitters each comprises an electrically conductive coil connected 
to a source of an alternating electrical current. 

A method according to claim 8 wherein the one or more magnetic 
emitters comprise at least three magnetic emitters. 

A method according to claim 9 comprising operating each of the at 
least three magnetic emitters to generate a corresponding magnetic 
signal having a frequency distinct from frequencies of other ones 
of the magnetic emitters. 

A method according to claim 9 wherein detecting the magnetic 
signals at six or more magnetic detectors comprises, for each of the 
at least three magnetic emitters, detecting the corresponding 
magnetic signal at each of a first set of magnetic detectors and the 
method comprises estimating a location of the magnetic emitter 
based upon relative amplitudes of the corresponding magnetic 
signal received at a subset of the first set of magnetic detectors 
which is different from a subset of detectors used for other ones of 
the at least three magnetic emitters. 
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12. A method according to claim 1 wherein determining relative 
amplitudes comprises using one of the magnetic detectors as a 
reference phase. 



13. A method according to claim 12 wherein using one of the 

magnetic detectors as a reference phase comprises computing: 



R; = Re 



i mO 



or a mathematical equivalent thereof wherein R q m is the relative 
amplitude of the signal from an q th one of the magnetic emitters 
10 detected by a m th one of the magnetic detectors; A q m is the complex 

amplitude of the signal from an q th one of the magnetic emitters 
detected by a m th one of the magnetic detectors; and A q m0 is the 
complex amplitude of the signal from an q th one of the magnetic 
emitters detected by a reference one of the magnetic detectors. 

15 

14. A method according to claim 1 wherein detecting the magnetic 
signals at six or more magnetic detectors comprises detecting the 
magnetic signals at each of a first set of magnetic detectors and the 
method comprises estimating a location of one of the magnetic 
20 emitters based upon relative amplitudes of the corresponding 

magnetic signal received at a subset of the first set of magnetic 
detectors. 



15. A method according to claim 14 wherein the first set of magnetic 
25 detectors includes at least 60 magnetic detectors and the subset of 
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the first set of magnetic detectors includes not more than 50 
magnetic detectors. 

16. A method according to claim 14 wherein the first set of magnetic 

5 detectors includes at least 150 magnetic detectors and the subset of 

the first set of magnetic detectors includes not more than 50 
magnetic detectors. 

17. A method according to claim 14 comprising selecting at a first 

10 group of magnetic detectors for inclusion in the subset based on an 

indicator of a variation in the magnetic signal with position of the 
magnetic emitter. 



18. A method according to claim 17 comprising selecting a second 
1 5 group of magnetic detectors for inclusion in the subset based upon 

a relative amplitude of the magnetic signal detected at the other 
magnetic detectors. 



19. A method according to claim 14 comprising selecting a first group 
20 of the magnetic detectors for inclusion in the subset based on a 

first criterion and selecting a second group of the magnetic 
detectors for inclusion in the subset based on a second criterion 
distinct from the first criterion. 



25 20. 



A method according to claim 1 9 wherein the first criterion selects 
for largest relative amplitude of the magnetic signal. 
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21. A method according to claim 20 wherein the second group selects 
for large variation in the magnetic signal with position of the 
magnetic emitter. 

5 22. A method according to claim 21 wherein the first group is larger 
than the second group. 

23. A method according to claim 14 wherein the subset includes 16 or 
fewer magnetic detectors. 

10 

24. A method according to claim 21 wherein the second criterion is 
based in part upon proximity to a reference magnetic detector. 

25. A method according to claim 21 wherein the second criterion 

1 5 comprises selecting magnetic detectors having relatively large 

values of V given by: 



V = 



\s — s 

I m max | 



r — r 

m max 



or a substantial mathematical equivalent thereof where s max is the 
amplitude of the signal at the magnetic detector having the largest 
20 signal from the magnetic emitter and f max is the location of the 

magnetic detector having the largest signal from the magnetic 
emitter. 



26. 
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A method according to claim 7 wherein the one or more magnetic 
emitters comprise at least two magnetic emitters for generating at 



-26- 

least two magnetic signals having frequencies f q which deviate 
from frequencies f qi at which the power spectrum of the detected 
noise signal has zeroes by at least two deviation amounts, and 
wherein performing the correction comprises calculating at least 
two corrected amplitudes, each of the at least two corrected 
amplitudes comprising a sum of a correction function applied to 
the at least two relative amplitudes. 

A method according to claim 1 wherein the one or more magnetic 
emitters comprise at least two magnetic emitters for generating at 
least two magnetic signals having frequencies f q which deviate 
from frequencies f qi at which the power spectrum of the detected 
noise signal has zeroes by at least two deviation amounts, and the 
method comprises correcting at least two relative amplitudes of the 
magnetic signals. 

A method according to claim 27 wherein correcting the relative 
amplitudes of the magnetic signals comprises calculating at least 
two corrected amplitudes, each of the at least two corrected 
amplitudes comprising a sum of a correction function applied to 
the at least two relative amplitudes. 

A method according to claim 28 wherein the correction function 
comprises a function of f q and f qg . 
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30. A method according to claim 29 wherein the correction function 
comprises a real correction function and an imaginary correction 
function. 

5 31. A method according to claim 30 wherein the interval extends from 
-T/2 to 772 and the real correction function comprises: 

2 T 

G c {f q ,f qi ) =-^ T cos(27f q t)e-' 2 ^'dt 

1 2 

sm(^f g -f qi )T) + sm[x(f g + f qi )T) 
\f q -f q )T + *{f q + f q )T 



and the imaginary correction function comprises: 
G s (f q ,f qi ) = - y j\sin(2^)e-' 2 ^ 

1 2 

"(fi-f*) T >{f*+f,)T 



10 32. A method according to claim 3 1 wherein the at least two magnetic 
emitters comprise three magnetic emitters for generating 
frequencies f l9 f 2 and f 3 which respectively deviate from 
frequencies f lh f 2i and f 3i at which the power spectrum of the 
detected noise signal has zeroes by first, second and third 

1 5 deviation amounts, and wherein calculating the corrected 

amplitudes from the relative amplitudes is accomplished by matrix 
equations: 
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G c (A,f u ) G c (f 2 ,f u ) G c (/ 3 ,/„) 
Get/:./,) G c (f 2 ,f 2i ) G c (/ 3 ,/ 2l ) 
Gc(/„/ 3 ,) G c (/ 2 ,/ 3 ,) G c (/ 3 ,/ 3/ ) 







X 


*2 







and 












G s (A,fu) 


G s (f 2 ,fu) 






G s {A,f 2i ) 


G s {f 2 ,f 2i ) 


A. 




G s (f„f») 


G s (f 2 ,f 3l ) 



where R } , R 2 , R 3 and 7 7 , 1 2 , I 3 respectively comprise real and 
5 imaginary parts of the relative amplitudes and C 7 , C 2 , C 3 and S h 

S 2 , S 3 respectively comprise real and imaginary parts of the 
corrected amplitudes. 

33. A method according to claim 28 where the integration time is an 
10 interval t 0 < t<t Q + T and the correction functions are: 



2 /o+7 ' 

G c {f q ,f qi )= ~ j 'cos(2tf q ty- n ^'dt = G^f q ,f qi )+G^f q ,f qi ) 

'o 

G s {f q >f gi )= f °)sm{27f q t)e- 2 *«'dt = G s R (f g ,f qi ) + G<(f q ,f gi ) 

'o 

and the matrix equation relating the real and imaginary amplitudes of the 
15 coil signals (R p and I p ) to the real and imaginary signals measured at the 
ideal frequencies (C p and S p ) is: 
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34. A method according to claim 28 wherein calculating the corrected 
amplitude comprises using a 3x3 matrix of complex numbers. 

5 35. A method according to claim 3 1 wherein calculating the corrected 
amplitudes comprises using a 3x3 matrix of complex numbers. 

36. A method according to claim 1 comprising determining positions of 
the magnetic emitters from the relative amplitudes of the magnetic 

10 signals. 

37. A method according to claim 36 comprising comparing the 
positions of the magnetic emitters to corresponding reference 
positions and triggering an alarm unless distances between the 

1 5 magnetic emitters and the corresponding reference positions are 

less than a threshold. 



38. A method according to claim 36 comprising comparing the 
positions of the magnetic emitters to corresponding reference 
20 positions and generating a sensory feedback indicator based upon 

departures of the positions of the magnetic emitters from the 
corresponding reference positions. 



